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Abstract 



Assuming the existence of one light sterile neutrino, we investigate the 
neutrino flavor mixing matrix in matter. Sum rules between the mixing pa- 
rameters in vacuum and their counterparts in matter are derived. By using 
these new sum rules, we obtain the simple but exact expressions of the ef- 
fective flavor mixing matrix in matter in terms of neutrino masses and the 
mixing parameters in vacuum. The rephasing invariants, sides of unitarity 
quadrangles and oscillation probabilities in matter are also achieved. Our 
model-independent results will be very helpful for analyzing flavor mixing 
and CP violation in the future long-baseline neutrino oscillation experiments. 
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I. INTRODUCTION 



The facts that neutrinos are massive and lepton flavors are mixed have been conflrmed 
by recent solar [1], atmospheric [2], reactor [3] and accelerator [4] neutrino oscillation exper- 
iments. However, the data coming from the Liquid Scintillation Neutrino Detector (LSND) 
experiment [5] have neither been confirmed nor refuted by the MiniBooNe experiment [6]. 
If the LSND result is proved to be true, most probably there exist more than three neu- 
trino mass eigenstates. From the phenomenological point of view, the most straightforward 
and attractive way to understand the result of LSND is to introduce the sterile neutrinos, 
which are electroweak singlets and do not couple to the Standard Model W and Z bosons. 
Such sterile neutrinos mix with the active neutrinos, thus they can significantly infiuence 
neutrino oscillations. The sterile neutrinos may also be the candidates of warm dark matter 
[8] in cosmology. Due to the Mikheyev, Smirnov and Wolfenstein (MSW) effect [9] within 
the supernovae, visible effects of sterile neutrinos in supernova explosion will probably be 
detected in the forthcoming neutrino telescopes [10]. There are also plenty of other reasons 
to believe the possible existence of sterile neutrinos [11]. Even if the LSND experiment is not 
confirmed, sterile neutrinos can not be fully excluded, and the number of neutrino species 
is still an important open question in neutrino physics. 

The direct effect of introducing sterile neutrinos is the mixing between active neutrinos 
and sterile neutrinos. For simplicity, here we only consider one light sterile neutrino. Hence 
the neutrino mass eigenstates will become linear combinations of three active neutrinos 
(i/g, 1/ , v^) and one sterile neutrino {Vg). To be explicit, we can write 
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where {i — 0, 1, 2, 3) denote the neutrino mass eigenstates, and V^j {a — s, e, /i, r) is the 
lepton fiavor mixing matrix. 

The behavior of neutrino propagation in vacuum can be described by the effective Hamil- 
tonian 

(VD'V^) , (2) 

where D = DiagjmQ, m^, mg}, m- are the masses of and E denotes the neutrino 
beam energy. In order to measure the lepton fiavor mixing angles precisely and to answer 
the question whether neutrinos violate CP, future accelerator experiments with very long 
baselines [12] are needed. In such long baseline experiments, neutrino beams should travel 
through the terrestrial matter, and their propagation can be significantly modified when 
they interact with matter. The effective Hamiltonian responsible for neutrinos travelling in 
matter can be written in a form analogous to Eq. (2) 

{VD'V^) , (3) 
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where D = Diag{?7io, ''Tii, '^sl; '"^j the effective neutrino masses in matter and V 
denotes the effective fepton mixing matrix in matter. Ti. is refated with Ti. as 

n-n = ^A, (4) 

where A — 2E x Diag(a', a, 0, 0) with the matter parameters a' — -\/2GpiV^/2 and a — 
\/2GpN^ arising from coherent forward neutrino scattering [13]. Here Gp is the Fermi 
constant, A^g and are the numbers of electrons and neutrons per unit volume. For most 
realistic long-baseline neutrino experiments, the background density can usually be treated 
as a constant [14] {N^ = constant and iV^ = constant). In this work, we take A^'g and A^^ as 
constants. 

In respect that the phenomenology of neutrino oscillations in matter can be formulated 

in terms of V and rh^, a number of authors have discussed the analytical relations between 
V and V [15-22] in the framework of three active neutrinos. The main purpose of this paper 
is to derive the exact analytical relationships between V and V in the four-neutrino mixing 
scheme. After an extension of our previous results [21] from three active neutrinos to three 
active neutrinos and one sterile neutrino, the model-independent sum rules between (V, rhj) 
and (y.m.^) arc given in section II. The relations between moduli \V^i\ and |V„j| are then 
obtained in section III, and so are the relations between Vg^jVg* and V^jV^j. Applications of 
our relations are also shown in section III. Finally, a brief summary is presented in section 
IV. 



II. SUM RULES IN THE FOUR-NEUTRINO MIXING SCHEME 

The effective Hamiltonians H and H in Eqs. (2) and (3) can be expressed as 

where the Greek indices a, 13 and so on generally run over s, e, and r, and the Latin indices 
i, j, etc. run over 0, 1, 2 and 3 here after. Eq. (4) turns out to be 

E ^ly^y^i = E mlV^i^;, + (6) 

i i 

with 

Here only the (s, s) and (e, e) elements of A are not vanishing. The unitarity conditions of 
V and V yield 
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These two sets of sum rules arc straightforward and obvious, but they are not enough for 
our calculations. In order to obtain the relations between V^jV^j and V^jVgj, we need more 
linear independent equations. 

Taking into account the unitarity of V{V), we find that = {1/2E)''VD^"V^ and 

(TY)" = (l/2E)"'i/L>2"yt hold, in which n = 0,±1,±2, etc. To be explicit, we obtain the 
elements of 'H{'H) as 

map ^ j^Y:m^''v^^v^, , 

map = j^i:rhrV,,V^i- (9) 

Note that the case n = corresponds to the unitarity conditions, and n = 1 corresponds to 
the sum rules in Eq. (6). For the case n = 2, we square the two sides of Eq. (6) and obtain 

i i 

In the calculation of Eq. (10), the expression of A in Eq. (7) has been considered. For the 
case n = 3, we get the forth set of equations by cubing both two sides of Eq. (6) 

i i 

Eqs. (6), (8), (10) and (11) constitute a full set of linear equations of V^^Vp^. By solving 
these equations, one can derive the exact expressions of V^jV^j in terms of mj, and V^jVg* . 
Note that, the sum rules obtained above are only valid for neutrinos. As for antineutrinos, 
the corresponding sum rules can directly be obtained by replacing [V , V, A) with [V* , V*, 
—A). Our formulae can be generalized to include more than one sterile neutrinos. In this 
paper, we only concentrate on the four-neutrino mixing scheme and the main results will be 
given in the next section. 



III. THE MIXING MATRIX ELEMENTS IN MATTER 

The sum rules obtained in the previous section totally contain ten sets of linear equations, 
which are related to the ten independent elements in 7Y. Let us first consider the case a — (3, 
in which the sum rules are simplified to 

y\v .p = y IV f = 1 , 

E^^lK.r = E^.'l^a.r + ^aa, 

i i 



4 



2 



\V -I 



+ "V rr?^!/ V ■V*V* A + 

' / > « 7j 71 aj "^77 1 "^aa " 



7>«J 



After a lengthy calculation, the solutions of Eq. (12) are obtained as 
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(13) 



(14) 



Here we have defined the mass-squared differences A^^ = mf — m^, A^,- = mf — m\ and 
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Aj = mf — vr?y In deriving Eq. (13), we have taken into account the relationship [23] 



(15) 



which can be clearly seen from the traces of two sides of Eq. (4). Note that only the 
mass-squared differences appear in our results, and the terms in proportion to the absolute 
neutrino masses disappear. This is consistent with the fact that oscillation processes are 
irrelevant to the absolute neutrino masses. Eq. (13) is very instructive and compact since 
it clearly shows the connections between and The survival probabilities denoted 

by P{ya ^a) or P{i^a ~^ ^a) [23] are only sensitive to the moduli of ViV) and thus 
these expressions will be very helpful for phenomenological analyses of the disappearance 
processes in the neutrino oscillation experiments. 

For the case a ^ (3, we obtain the explicit expressions 
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and a ^ /3. In the hmit ^ — > 0, our results are reduced to the vacuum case V ^V. 

To show the applications of our results, we define the rephasing invariants in matter as 
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(18) 



which are similar to the one in the three-neutrino mixing scheme [24]. Only nine of these 
invariants are independent [25]. We can choose any set of nine independent J^-'^ to fully 
describe CP and T violation in matter. 

The twelve unitarity quadrangles [25], which are the geometrical presentations of the 
orthogonality of V{y) in the complex plane, are well related with CP and T violation in 
neutrino oscillations. The twenty-four sides of six quadrangles have been given in Eq. (16). 
Once precise measurements are provided in future, we can systematically analyze the matter 
effects on these quadrangle sides by taking use of our relations. 

The probability of one neutrino converting to another neutrino in matter is given 

as 

P{^a - = <5a/5 - 4 E {v^^Vp^V:^V;) siu' F,,] - 2^^ siu 2F,,) , (19) 

i<j i<j 

where Fji = 1.27 A j^L/ E with L standing for the baseline length (in unit of km) and E 
being the neutrino beam energy (in unit of GeV). Prom our simple and exact formulae, the 
oscillation probabilities for all channels in matter can be calculated and the CP-violating 
information can also be extracted. 

It should be mentioned that our formulae are also very suitable for numerical analyses, 
and the relevant work will be elaborated elsewhere. 



IV. SUMMARY 

Motivated by the interesting extension of the standard flavor mixing scheme of three 
active neutrinos with one sterile neutrino, we have derived a full set of sum rules in the four- 
neutrino mixing scheme. By using these sum rules, we have calculated the effective neutrino 
mixing matrix in matter V in terms of m^, m- and the genuine neutrino mixing matrix V. 
The rephasing invariants, sides of six unitarity quadrangles and oscillation probabilities in 
matter can be obtained by using our formulae in Eqs. (13) and (16). These compact and 
simple results are expected to be very helpful to describe the leptonic flavor mixing and CP 
and T violation in matter in a model-independent way. The relations are also very useful 
for a systematic study of neutrino oscillations in the future long-baseline experiments. 

The author would like to thank Professor Z.Z. Xing for reading the manuscript, making 
many corrections and giving a number of helpful suggestions. He is also grateful to S. Zhou 
and W.L. Guo for useful discussions. This work was supported in part by the National 
Nature Science Foundation of China. 
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